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Abstract 

Cosmological perturbation equations derived from low-energy effective ac- 
tions are shown to be invariant under a duality transformation reminiscent 
of electric-magnetic, strong- weak coupling, S-duality. A manifestly duality- 
invariant approximation for perturbations far outside the horizon is intro- 
duced, and it is argued to be useful even during a high curvature epoch. 
Duality manifests itself through a remnant symmetry acting on the classical 
moduli of cosmological models, and implying lower bounds on the number 
and energy density of produced particles. 
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I. INTRODUCTION 



In inflationary cosmology JT| the Universe quickly becomes homogeneous, isotropic, and 
spatially flat, but quantum mechanics superimposes on the smooth classical backgrounds 
calculable inhomogeneity perturbations [0. It is the purpose of this paper to point out 
a general duality symmetry obeyed by the lowest-order equations for these perturbations, 
and to show that duality may be used to obtain lower bounds on the energy density of the 
produced particles. 

Our analysis applies to any physical perturbation \1/ that is described, to lowest order in 
the strength of the perturbation and in its derivatives, by the following quadratic action 



drj d A x S(rj) m' 2 - (Vtf ) 2 . (1.1) 



In eq. ( |1.1|) , which is generic in low-energy spatially flat backgrounds, rj is the conformal time 
coordinate, and a prime denotes d/dr/. The function S(rj), sometimes called the "pump" 
field, is, for any given a known function of the scale factor, a(t]), and of additional scalar 
fields (such as a dilaton 4>(r]), moduli bi(rj), etc.) on which the background may depend. 

The duality symmetry, which, to the best of our knowledge, has not been noticed before in 
the context of cosmological perturbation theory, involves the interchange of the perturbation 
\l/ with its conjugate momentum II, as well as the inversion of the pump field S. In some 
cases it corresponds to a time-dependent version of T-duality, interchanging small and large 
scale factors HH; in others it corresponds to a time-dependent version of ^-duality [HH, 
interchanging weak and strong coupling and electric and magnetic degrees of freedom; in 
most cases it amounts to a time-dependent mixture of both. As in other applications, duality 
turns out to be a powerful tool for obtaining results that are hard to obtain otherwise. 

Although most of our discussion will be completely general, it may be useful to recall that, 
in the context of string cosmology, the massless perturbation \1/ may represent a physical 
polarization of tensor metric perturbations (i.e. a gravity wave), physical electromagnetic 
perturbations (i.e. photons), axionic perturbations, and so on. In these three cases the pump 



2 



field S corresponds, respectively, to the Einstein-frame scale factor given in terms of the 
dilaton as S = a 2 e~^, to the inverse string coupling || S = e~^, and to the "axion-frame" 
scale factor S = a 2 e^. Here, and in what follows, a is the scale factor of the so-called 
string-frame metric, which coincides with the usual scale factor of general relativity whenever 
the dilaton is constant. For a generic string cosmology perturbation S(rj) = a 2m e l ^b 2m \ 
where the numbers /, m, rrii depend on the type of perturbation. Arguably, the action ( |1 . 1|) 
provides a useful description for a much wider class of perturbations. 

II. DUALITY 

In order to discuss the duality symmetry of the action (|1 . 1|) , it is convenient to use the 
Hamiltonian formalism. The Hamiltonian corresponding to ( |1 . 1|) is given by 



H = - I <r\r 



{2.1) 



where II is the momentum conjugate to The first-order Hamilton equations read 



= ^ = S-'U , U' = ~ = SV 2 ^ , (2.2) 

and lead to the decoupled second-order equations 

^" + ^'-v 2 ^ = o, n"- ^n'- v 2 n = o, (2.3) 

of which only the first is commonly used in the study of perturbations. In Fourier space the 
Hamiltonian ( |2.1|) is given by 



(2.4) 

where = and n_£ = Tit. The equations of motion become 

% = S-'K-k , = -Sk^_j; , (2.5) 

where k — \k\. The transformation 
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ns-n s =** s> ^s-*^-*"^, 5- SW" 1 , (2.6) 

leaves the Hamiltonian, Poisson brackets, and equations of motion unchanged. 

We shall call this transformation, for short, S'-duality since it contains the usual strong- 
weak coupling duality as a special case. Indeed, by perturbing the tree-level, lowest- 
derivative string effective action, one finds that, for heterotic, four-dimensional gauge bosons, 
the function S is given simply by |J S(rj) = e~^ v \ Since determines the effective string 
coupling constant, = g 2 string1 the transformation S — > S = S -1 amounts to performing, 
at each time 77, a strong- weak coupling transformation g string — > g^trmgi which is part of the 
SL(2,Z) S-duality group ||. Moreover, in the case of electromagnetic perturbations, the 
contribution of to the Hamiltonian can be identified with that of the magnetic field, 
whereas IT is proportional to the electric field, hence, in this case, the transformation ( |2.6|) 
exactly coincides with electric-magnetic duality. 

We are interested in a situation where the time evolution of the pump field amplifies the 
zero-point vacuum fluctuations of the field satisfying the initial conditions 

(S^II 2 ) = (S(W) 2 >, (2.7) 

where (• • •) denotes ensemble average, or expectation value if perturbations are quantized. 
Under S-duality the equations of motion are form-invariant, while the initial conditions are 
invariant on the average in the sense of eq. ( |2.7|) . Therefore, under S'-duality, solutions of 
the perturbation equations properly normalized to zero-point fluctuations, are mapped into 
other normalized solutions with the same total Hamiltonian. We will show that, as a con- 
sequence of this property, the energy spectrum can be consistently estimated by truncating 
the solutions outside the horizon to the frozen modes of \I/ and n. 

This may be regarded as trivial, if one has in mind a background where the non-frozen 
part of \I/ quickly decays outside the horizon, as in the standard inflationary scenario [0. 
There are backgrounds, however, in which \l/ grows in time instead of decaying, and the 
growing mode of \l/ dominates the perturbation spectrum. Typical examples can be found 
in scalar-tensor models of gravity and, in particular, in string cosmology DTTJ. In this context, 



the possibility of a consistent truncation to the frozen part of the perturbations is a powerful 
result. Only in a duality- invariant approach is such a truncation consistent, since, as shown 
below, the contribution to the Hamiltonian of the growing mode of $ is simply replaced by 
the contribution of the frozen part of its duality-related conjugate momentum. 

III. APPROXIMATE SOLUTIONS INSIDE AND OUTSIDE THE HORIZON 

In order to solve the perturbation equations, and to normalize the spectrum, it is con- 
venient to introduce the canonical variables II, whose Fourier modes are defined by 

^ k = s^ 2 ^ k , n fc = s~ 1 ' 2 u k , (3.1) 

so that the Hamiltonian density takes the canonical form 

H= 1 -Y,{\n k \ 2 + k 2 \$ k \ 2 ). (3.2) 

k 

Under S'-duality these canonical variables transform as the original variables. They satisfy 
Schrodinger-like equations: 

$ k " + [k 2 - (s^y's- 1 ? 2 ] § fc = o , n fc " + [k 2 - (s-^ys 1 / 2 ] n fc = o. (3.3) 

The amplification of perturbations is typically associated with a transition from an infla- 
tionary phase in which the pump field is accelerated (in cosmic time), to a post-inflationary 
phase in which the pump field is decelerated or constant. In such a class of backgrounds, 
the "effective potentials" V* = (S 1/2 )" 'S~ 1 / 2 , Vfi = (S- 1 ' 2 )" 'S 1 ' 2 , grow (in absolute value) 
during the phase of accelerated evolution, and decrease in the post-inflationary, decelerated 
epoch, vanishing asymptotically for very early times 77 — > —00 and vanishing again for very 
late times 77 — > +00. 

The initial evolution of perturbations, for all modes with k 2 > \Vq,\, \Vn\, may be de- 
scribed by the WKB-like approximate solutions of eqs. (|3.3|) 



- , , , 2 v-l/4 -*J>'(* 2 -^) V2 

*k(v) = (k 2 -V*) e •» 



,1/2 



, 9 N -i/4 -JM* 2 -^) 

U k (r ] ) = k(k 2 -V u ) e "0 , (3.4) 



which we have normalized to a vacuum fluctuation, and where the extra factor of k in the 
solution for 11^. comes from consistency with the first-order equations ( |2.5| ). We have ignored 
a possible relative phase in the solutions. Solutions ( |3.4| ) manifestly preserve the S-duality 
symmetry of the equations, since the potentials V®, Vn get interchanged under S — > S -1 . 

In the case of power- law evolution of the pump field, S 1 ^ 2 ~ r/ a , the effective potentials 
grow like r]~ 2 as 77 — > 0, and the approximate solutions (|3.4|) are valid, at a given time 77, 
for all modes k with (krf) 2 > 1. These modes are usually referred to as being "inside the 
horizon". The corrections to solutions (13.41) are of order V^j(k 2 - V*), V*"/(k 2 - Vy) 
and higher, which are indeed small for large (kt]) 2 . Asymptotically, at 77 — > —00, we have 



— > and solutions Q3.4Q reduce to correctly normalized vacuum fluctuations 

$ h (ri) = k~ 1/2 e , n fc (77) = k +l ' 2 e ~ ikr >+% : (3.5) 

where ip^, ip'^ are random phases, originating from the random initial conditions. Note that 
because of the random phases, ^-duality holds only on the average in the sense of eq. ( j2.7|) . 

For k 2 < \V\p\, \Vu\, equivalently [kr]) 2 < 1 if S 1 / 2 ~ 77", the modes are usually referred 
to as being "outside the horizon" . For such modes, it is possible to write "exact" solutions 
to eq. ( p.3| ) as follows, 

$ fc (7j) = S 1/2 A k Tcos(S-\ S) + B k Tsin(5~ 1 , S) 
n k ( v ) = kS- 1/2 B k Tcos(5, S- 1 ) - A k Tsin(S, S' 1 ) , (3.6) 

where A k ,B k are arbitrary integration constants, and the functions Tcos(5~ 1 , S; k, 77), 

Tsin(S' _1 , S; fc, 77) are defined through the following expansions: 

m 



Tcos(S 1 , S; k, 77) = 1 — k J dr]iS 1 (tji) k J drjzSfa 

T) r)i r/2 ?73 

+ k J c/77 1 S ,_1 (t7i) k J dr]2S(r] 2 ) k J drjsS' 1 ^) k J di] 4 S(r] 4 ) H ; 

Vex ^)ex Vex "tyex 

V V Vi 12 

Tsin(S' _1 , 5"; k,rj) = k J drjxS" 1 (771) - k J dr^S' 1 ^) k J drftSfa) k J drjsS" 1 (r} 3 ) 

Tjex "Hex f)ex 1]ex 

+ ■••■ (3.7) 



The time r/ ex , appearing as an arbitrary lower limit of integration in ( |3.7|) , is most conve- 
niently chosen to be near horizon-exit time, i.e. such that krj ex ~ 1. The above functions 
satisfy 



Tcos(S-\S) 



Tsm(S-\ 5)1' = | Tcos(S, S' 1 ). (3.8) 



Using these relations, and similar ones for [Tcos(5', S" -1 )]' and [Tsin(S l , 5 1 " 1 )]', it is possible 
to verify explicitly that ( WM are indeed solutions of eqs. ( |3.3|) . Formally, these solutions 
are also valid inside the horizon, but the power series defining Tcos, Tsin are not convergent 
there. Solutions (|3.6| ) manifestly preserve S'-duality, which now acts on the parameters of 
the solutions in the following simple way: 



-B k 



Bk — > A k 



S -> s- 1 . 



(3.9) 



We now fix the integration constants by matching solutions ( |3.4| ) and (|3.6|) and their first 
derivatives at horizon-exit time, for which kr] ex ~ 1. The result (written for the variables 
without the hats) is simply 



V k ( V ) = (k 2 - V£ 
Mv) = k(k 2 - Vg 



-1/4 



-1/4 



S~ 1/2 Tcos(S-\ S)e 1 ^ + S l J 2 Tsm{S~\ S) e i<p S 
SH 2 Tcos(S, S- 1 ^ - S; x 1/2 Tsin(S, S^e*** 



(3.10) 



where tp^ and tp'* now contain an additional phase depending on the exit time i] ex , and a 
label ex means evaluation at rj = r] ex . The above matching respects S'-duality on the average, 
in the sense of eq. 



IV. ENERGY SPECTRA 

To complete a calculation of the energy spectrum of the amplified perturbations we need 
to perform a second matching step, to late-time solutions after horizon re-entry in the decel- 
erated, post-inflationary era. However, if we are only interested in the leading contribution 
to the final energy spectrum, we can estimate it by using the well-known phenomenon of 
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freezing of perturbations. Solutions ( |3.10|) always contain a time-independent part coming 
from the first term of the Tcos expansion ( |3.7| ), and always either the constant part of \& or 
the constant part of II dominates the Hamiltonian in the large wavelength limit \krj\ <C 1. 

In order to illustrate this important point, we parametrize the background during the 
phase of accelerated evolution with a power-like behaviour of the pumping field, S 1 ^ 2 ~ \r/\ a 
for i] < 0, 7] — > 0_. From the general solution (|3.6|) we obtain, in the limit \kr)\ -C 1, 

$ k = A k \v\ a + ^-B k \v\~ a , fi fc = k (B k \rj\- a - J^LA k \ri\ a j , (4.1) 
1 — 2a \ 1 + 2a I 



with logarithmic corrections for a = ±1/2. Computing the Hamiltonian density ( p.2| ) we 
obtain, for \kr]\ <C 1, 

H = \ Y. k4\A k \ 2 \r,\ 2 ° [l + 0(|*r,| 2 )] + |B t | ! M _2 ° [l + 0(|Ar,| 2 )] + O (\A k \\B k \\k n \) 1 



k K 



m v\ 2a + \B k \ 2 \v\~ 2a \- (4-2) 



Thus, irrespectively of the value and sign of a, the joint contribution of the constant modes 
~ A k and 11^ ~ kB k always provides an accurate estimate of the leading contribution to 
the Hamiltonian for super-horizon wavelengths. Which of the two constant modes in ( [1.2| ) 
gives the dominant contribution varies from case to case. Also, it is not always true that 
the subleading constant mode provides the leading correction. 

The above statements can be explicitly checked when the pump field 5* is given by a power 
of 7] with different powers in the different cosmological epochs. In this case, the joining of 
solutions at the transition times can be done analytically, exact formulae for the Bogoliubov 



coefficients can be given, and their duality can be directly checked [12]. An important 



point Jl2| is that, for backgrounds with discontinuous derivatives at the transitions, exact 



duality is satisfied only if continuity is imposed on the original perturbation (and on its first 
derivative) and not on the rescaled fields n. 

Up to numerical factors, the contributions of V^ x , V^ x in ( |3.10| ) can be neglected, and one 



can thus approximate the normalized vacuum fluctuations outside the horizon by truncating 
to 1 the Tcos expansion and to the Tsin expansion, 
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Vk(v) = (kS ex r 1/2 e^ 



n k (v) = (ks, 



1/2 up' 
ex J e K ■ 



(4.3) 



This truncation reflects just a remnant of S-duality, the transformation S ex — > S~ x . It is 
adequate for an approximate computation of the energy spectrum after re-entry. 

At late times 77 — > +00, when V$, Vu — > in the post-inflationary phase, the perturbation 
modes with ^ > 1 oscillate freely inside the horizon, according to the equations 



§ fc " + k 2 $ k = , 



fl fe " + k 2 YL k = 0, 



(4.4) 



with plane-wave solutions 



n fc = % 



2S 



iVfc (-C k e- ihn + D*e 



(4.5) 



The matching of these solutions to the approximate solutions ( [4.3| ), performed at 77 = r\ T 
k~ l , imposes the conditions 



1 / S PX \ / 



e *Pk 



Q x 1/2 , 



VP-, 

e k 



(4.6) 



from which we obtain C k , D k . The final amplified vacuum fluctuation amplitude, in the 
regime krj ^> 1, is thus given by 



g x -1/2 . / s ex \ 1/2 i , 

— ^ J cos(krj) e lip * + I — ^ j sin(/c^) e i</3 * 

Sex\ 1/2 i ' ( S ex \~ 1/2 

— j^- ) cos(kr]) e tLp k — I — J suffer/) e lip * 



(4.7) 



(we have absorbed into (p k and <^' fc an additional phase factor kr] re , arising from the matching). 
The corresponding averaged Hamiltonian density, for a mode k, 



(^) = ^((|n fc | 2 ) + ^ 2 (l^ 



, / S ex S re 
M^- + ^-|' 



(4.8) 



is invariant under S ex — > S~ x l , S re — > S^, and also under an overall rescaling of 5*. The 
same invariances thus characterize the spectral energy distribution, dpk/dln k = (k 3 /a 4 ){H k ) , 
which can be written in terms of the proper energy u = k/a as 



dp(u) _ 4 



din a; 



S ex (uj) S re (uj) 



S re (uj) S ex (u) 



c/Max{|p,§4. (4.9) 

l_ >-Ve O ex j 



The approximate expression ( }4.9|) reproduces, up to numerical factors 0(1), the leading 
contribution to all known particle spectra for various models of cosmological evolution. 
Consider, for instance, graviton production in a transition at rj = 77* from de Sitter inflation 
to radiation, at constant dilaton. In that case S = a 2 , S ex / S* = (krj^) 2 = (k/k*) 2 and 
S*/S re = (fc/fc*) 2 , where k* = 777 1 is the maximal amplified frequency. Therefore S re /S ex = 
(fc*/fc) 4 . For all modes u < a;*, where oj* = k*/a, M&x{S ex / S re , S re / S ex } = S re /S ex = 
(u^/cu) 4 , and we recover the well-known flat Harrison- Zeldovich spectrum ||, dp/dkiui ~ cu 4 . 

If, instead, we consider the transition from a phase of dilaton-driven inflation to radi- 
ation, typical of string cosmology, we can recover from (|4.9|) the leading contribution to 
known spectra for photons 0, gravitons and axions [|T^,[y||l2| . For heterotic photons, in 
particular, and more generally for gauge bosons coming from extra dimensions, S(r]) = e~^ v ^ 
coincides with the real part of the usual S'-field of string theory, and the leading contribu- 
tion to the energy spectrum ( |4.9| ) is the same for two cosmologies that arrive at the same 
final S from above (weak coupling) or from below (strong coupling). In the two cases the 
energy is predominantly stored in the electric and in the magnetic field, respectively. From 
the approximate expression ( |4.9|) we can also recover the leading contribution to the energy 



spectrum for more general classes of backgrounds and types of perturbations [13.12 



V. LOWER BOUNDS ON ENERGY SPECTRA 



Since the approximate expression for the energy spectra, eq. ( f4.9| ), contains, as a direct 
result of the remnant duality symmetry, a factor of the form x + an interesting conse- 
quence of the remnant duality symmetry is the existence of a uniform and model-independent 
lower bound for the energy spectrum: 

* ^ (5-1) 

a In a; 
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valid for all models of backgrounds and all types of perturbations: it represents the minimal 
spectrum corresponding to the production of one particle per mode, per polarization and 
per unit phase-space volume ||14|| . For some particular model of background, and for a given 
type of perturbation, eq. ( |4.9| ) can provide a less trivial lower bound, corresponding to a 
"minimal" spectrum flatter than u 4 . In such cases the frequency dependence of the lower 
bound is, in general, model-dependent. 

Consider, in particular, the amplification of tensor perturbations in a string cosmology 
background in which the initial dilaton-driven phase, with 5* ~ rj, is followed at rj = r] s by a 
high-curvature string phase, with final transition to the radiation-dominated era at r] — rji. 
During the string phase, a total redshift z s = \a\/a s \ > 1 is accumulated, the curvature is 
approximately constant, and the dilaton grows by the ratio exp(0! — 4> s )/2 = g\/g s > 1. 
Equation (fO|) gives, in that case: 

dp -^Maxl^^V^-f^V^ 3 ), uxu.^^/z.. (5.2) 



dhauj { uj \gi J s ' oji \g. 

The assumption that the pump field keeps growing during the string phase, z s g s jg\ > 1 
then leads to the well-known cubic (n = 4) graviton spectrum M, whose lower bound, 



dp ^ 4 (9s\ 2 ( u 



3 



7 1 £ ^ - - , (5-3) 

is more constraining than eq. (|5.1|) . The "minimal" spectrum is dp/d\nuj = UiU 3 and 
corresponds to the limiting case in which g s — > gi and z s — > 1, namely to a model with almost 
no intermediate string phase. Note that, for gravitons, the largest sub leading contribution is 



not given by the second term in the bracket in eq. ( |5.2|) |8|JT3|] but, instead, can be extracted 
by carefully matching solutions at rj — r\\. When this is done, the more stringent bound can 
be seen as a consequence of a remnant symmetry acting on both z s and g s /gi- 

An important point is that the duality-symmetric form ( |4.9| ) of the spectrum, which 
provides an accurate order-of-magnitude estimate for the energy spectra produced from per- 
turbations whose action is given by ( |1 . 1| ) , can also provide a lower bound on the total energy 
density of long-wavelength perturbations, even when the background evolution includes a 
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high-curvature phase in which the action and perturbation equations are not known. The 
argument goes as follows. 



The result ( [4.9| ) was obtained by approximating the evolution of \I> and II by their 
constant modes outside the horizon. The approximation was motivated by the explicit 
solution of the perturbation equations. When the background enters a high- curvature phase, 
which requires a corrected effective action for an accurate description, the exact form of the 
perturbation equations and their solutions are not always known, or in any case difficult to 
obtain in practice. We expect that, for any type of action, the perturbation equations should 
always admit constant solutions whenever spatial gradients are negligible with respect to 
time derivatives, in agreement with the general physical principle of freezing of perturbations 
outside the horizon. 

For metric perturbations, a formal argument can be given to justify our conjecture. Since 
the action is invariant under gauge (i.e. coordinate) transformations, a particular set of al- 
lowed perturbations are infinitesimal gauge transforms of the background itself. Of course, 
these are spurious perturbations without any physical content. It is possible, however, to 
construct physical, large-wavelength perturbations, which approach the infinitesimal gauge 
transformations in the infinite- wavelength limit. By continuity, these large- wavelength grav- 
itational waves are approximate solutions of the perturbation equations. The amount by 
which they deviate from exact solutions away from the limit is controlled by the only di- 
mensionless parameter present, i.e. by kr]. We may thus conclude that, at least for metric 
perturbations, the constant mode in \1/ is always present for perturbations sufficiently outside 
the horizon. Canonical Poisson brackets then require that the same be true for II. Although 
we cannot exclude that other modes in \I/ and II will dominate over the constant modes, 
assuming that this does not happen will provide a lower bound on the amplified spectra in 
the presence of a high-curvature phase. 

For tensor perturbations, this argument has been confirmed within an explicit string 
theory model of quadratic curvature corrections [|T5|]. The results of a numerical integration 
show that \1/ becomes frozen outside the horizon even after including the higher-curvature 
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terms in the perturbation equations. The final normalized amplitude is larger, by a numerical 
factor, than the estimate ( |4.3|) , so that the spectrum ( |4.9| ) indeed represents a lower bound 
for the total energy density. 

VI. CONCLUSION 

We have shown that the equations of cosmological perturbation theory, to lowest order in 
the strength of perturbations and in their derivatives, are invariant under a duality transfor- 
mation acting on the classical moduli of the background and interchanging the perturbation 
with its conjugate momentum. As a consequence of this invariance, the large-wavelength 
part of the Hamiltonian is always dominated by the sum of the contributions of the frozen 
modes of the perturbation and of its conjugate momentum. 

The energy spectrum obtained by truncating the solutions of the perturbation equa- 
tions to the constant modes is characterized by a residual duality symmetry; it reproduces 
known results of produced particle spectra and can provide a useful lower bound on parti- 
cle production, when our knowledge of the detailed dynamical history of the background is 
approximate or incomplete. 
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